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1.  Introduction.  In  many  situations  where  the  lifetimes  of  systems  and  components 

« ♦ 

are  considered,  these  lifetimes  do  not  have  exponential  distributions.  Rather, 

% 

these  lifetimes  reflect  the  effect  of  wearout.  Four  alternative  classes  of 
lifetimes  which  describe  various  types  of  wearout  and  have  been  extensively 
studied  are:  1)  the  increasing  failure  rate  (IFR)  class;  2)  the  increasing 
failure  rate  average  (IFRA)  class;  3)  the  new  better  than  used  (NBU)  class;  and 
4)  the  new  better  than  used  in  expectation  (NBUE)  class.  In  the  case  where  the 
component  lifetimes  can  be  assumed  to  be  independent  these  classes  have  proven 
to  be  very  useful.  Recently  in  an  attempt  to  describe  the  more  realistic 
situations  where  this  independence  assumption  cannot  be  made,  various  multi- 
variate versions  of  the  above  classes  have  been  proposed. 

In  Section  2 we  shall  discuss  the  univariate  classes  and  their  properties. 

Then  we  shall  discuss  the  most  important  multivariate  analogs  of  the  IFR  (Section 
3)  and  the  IFRA  (Section  4)  classes.  Properties  of  these  classes  as  well  as 
their  relation  to  other  proposed  classes  will  be  given.  In  Section  5,  new 
multivariate  definitions  of  NBU  and  NBUE  classes  will  be  given  and  these  will  be 
discussed. 

Throughout  the  paper,  the  terminology  and  notation  of  Barlow  and  Proschan 
[3]  will  be  used  with  the  exception  that  a structure  function  <J>(x)  will  be 
called  monotone  if  it  is  nondecreasing  in  each  of  its  components  and  in  addition 
$(0)  = 0 and  <p  (1)  = 1. 

2.  Univariate  Classes.  The  most  important  and  most  studied  class  of  lifetimes 
which  describe  wearout  is  the  class  of  distributions  with  increasing  failure 
rate.  This  concept  has  been  used  in  actuarial  science,  in  statistics,  and  in 
engineering  reliability  where  it  is  sometimes  called  "increasing  hazard  rate". 

Let  T be  a random  lifetime  with  distribution  function  F(x)  * P{T  < x}  and 


F(0)  - 0*.  survival  function  F(x)  - 1 - F(x)  and  density  f(x)  (if  it  exists). 
Then  T (or  F)  is  said  to  have  increasing  failure  rate  (IFR)  if 


P(x  < T s x + t) 
P{ x < T) 


increases  in  x > 0 for  t > 0. 


(2.1) 


This  means  that,  given  that  T has  survived  beyond  time  x,  the  probability  that 
it  fails  within  the  next  t units  of  time  increases  as  x increases,  i.e.  the  older 
the  lifetime,  the  more  likely  it  is  to  fail.  A lifetime  is  said  to  have 
decreasing  ,railure  rate  (DFR)  if  the  quantity  in  (2,1)  decreases  in  x ' 0 for 
t >_  0.  The  definition  of  IFR,  more  often  encountered,  which  is  an  instantaneous 
version  of  (2.1)  and  is  equivalent  to  it  if  the  density  exists,  is  that  the 
failure  rate  function 


r(x)  “ — — increases  in  x > 0.  (2.2) 

V(x) 

A corresponding  statement  holds  for  DFR.  Another  way  of  writing  (2.1),  which 
is  useful  in  the  multivariate  cases  to  be  described,  is  that 


F(x  + t)  _ P{T  > x + t) 
F(x)  PIT  > x) 


1- 


P(x  < T ^ x + t) 
P{x  < T) 


decreases  in  x > 0 for 


t > 0. 

(273) 


This  amounts  to  saying  that  the  survival  rate  decreases  as  the  age  of  the 
component  increases.  Similarly  an  increasing  version  of  (2.3)  (i.e.  increasing 
survival  rate)  is  equivalent  to  DFR.  Still  another  version  of  IFR  (DFR)  is  that 
the  function 

log  F(t)  is  concave  (convex).  (2.4) 

The  classes  of  IFR  and  DFR  lifetimes  satisfy  the  following  properties: 

* 

The  usual  assumption  is  that  F(0~)  * 0.  We  make  the  simplifying  assumption 
F(0)  - 0 here  for  the  purpose  of  this  exposition.  It  lends  to  slightly  simpler 


def inltions . 
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l)  the  only  lifetimes  which  are  both  IFR  and  Dl'R  are  the  exponential  lifetimes; 
li)  convolutions  of  TFR  lifetimes  are  IFR  (this  means  that  for  independent  IFR 
lifetimes,  if  a failed  component  is  replaced  by  a spare,  then  the  accumulated 
lifetime  is  IFR).  See  Barlow  and  Proschan  [3]  for  additional  properties  and 
discussion. 

The  second  major  class  of  lifetimes  we  will  discuss  is  the  class  whose 
lifetimes  have  increasing  failure  rate  average  (l FRA).  These  are  the  lifetimes 
whose  survival  function  satisfies 

F(at)  > F^t)  for  all  0 < a < l and  all  t > 0.  (2.5) 

If  the  density  exists  it  is  not  hard  to  show  that  this  is  equivalent  to  the 
more  intuitive  condition  that 

J*  r (x)dx 

increases  in  t > 0 (2.6) 

where  r(x)  is  the  failure  rate  function  of  (2.2),  and  so  (2.6)  gives  that  the 
average  failure  rate  increases.  Another  equivalent  form  of  1FRA,  which  was  used 
by  Block  and  Savits  *5]  to  establish  that  the  convolution  of  IFRA  lifetimes 
is  IFRA,  is  that 

/CO  on 

hn(x/n)  dF(x)  > {/_  h(x)  dF(x) l , all  0 < « < 1 (2.7) 

0 ~ '0  — 

and  all  nonnegative  nondecrcasing  functions  h. 

The  class  of  IFRA  lifetimes  which  is  wider  than  (i.e.  contains)  the  class 
of  IFR  lifetimes  satisfies  the  following  properties:  i)  for  monotone  systems 
with  independent  IFRA  lifetimes,  the  system  lifetime  is  IFRA;  ii)  the  IFRA  class 
is  the  smallest  class  containing  the  exponential  distributions  which  is  closed 
under  the  formation  of  coherent,  systems  and  limits  in  distribution  (see  Ross  [24] 
for  an  extension  to  IFRA  processes);  iii)  the  lifetime  of  a device  arising  as  a 
consequence  of  a natural  nonfatal  shock  model  is  IFRA  (see  Theorem  3.8,  p.  94, 


Barlow  .-mil  Broschan  |11);  iv)  convolutions  of  IFRA  lifetimes  are  1FRA  (see 


Mock  and  Savlts  [■>])•  For  other  properties  of  this  class  t no  hull  up,  statistical 
procedures  see  Barlow  and  Marshall  [21,  Barlow  |1],  and  Doksnm  [131,  (141. 

Still  wider  classes  of  lifetimes  are  those  which  are  new  better  than  used 
(NBU)  and  new  better  than  used  in  expectation  (NB1IK) . A lifetime  T is  NRU  if 
the  probabilitv  that  it  will  survive  when  new  is  greater  than  the  probability  that 
it  will  survive  given  that  it  is  anv  other  age,  i.e. 


x + t 


- -FXp-(V-)-tj  < F(t)  - 1MT  tl  for  all 


(?.R) 


A lifetime  is  NBl'K  if  an  integrated  version  of  (2.7)  is  satisfied,  i.e. 

/«»  — — 

F(x)  dx  ' u Fit)  for  t ' 0 (2.^) 

t 

.«*  - 

where  p ■ ) F(x)  dx  is  assumed  to  be  f inite. 

Another  wav  of  expressing  that  a lifetime  T is  NRIIF.  is  that 

K(T  - t | T ' t)  < H for  t ' 0 (2.10) 

which  means  that  given  anv  fixed  age  the  residual  mean  lifetime  is  smaller  than 
the  mean  lifetime. 

Both  the  NBU  and  the  NBUF.  classes  have  been  shown  to  be  useful  in  the 
solution  of  certain  maintenance  problems.  In  particular,  lifetimes  are  NBU  if 
replacement  policies  are  to  be  beneficial  in  a certain  sense.  The  NRUK  life- 
times are  the  largest  class  for  which  the  number  of  failures  observed  in  a 
"replace  at  failure  only"  policy  is  larger  stochastically  in  the  long  run  then 
when  the  process  starts.  See  Chapter  6 of  Rarlow  and  Proschnn  [3]  for  details. 
Another  situation  where  NBU  arises  is  given  bv  Ross  [231  who  shows  that  for  a 
monotone  system  with  independent,  exponential  components  and  exponential  repair 
times,  the  tins'  until  first  system  failure  is  NBU.  Furthermore,  monotone  systems 


i 

- 


- 
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with  independent  NBU  components  have  NBU  lifetimes  and  convolutions  of  NRU 
(NRUE)  lifetimes  are  NBU  (NBUE) . Tests  for  NBU  and  NRUE  have  heen  given  bv 
Hollander  and  Prosehan  [18],  [191. 

3.  Multivariate  IFR.  There  have  been  several  different  multivariate  extensions 
of  the  concept  of  increasing  failure  rate.  Most  of  these  have  heen  based  on  one 
of  two  notions.  The  first  of  these  arises  from  a generalization  of  (2.3)  and  a 
second  from  a generalization  of  (2.4). 

Let  T * (T  , ....  T ) be  a multivariate  lifetime  and  let  F(t)  = F(t, t ) 
la  —In 

- P{T  > t .....  T > t ) be  the  survival  function.  Also  assume  that 

11  a n 

F(0)  » 1 where  0 = (0,...,  0).  Since  (2.3)  is  interpreted  as  a device  being 
less  likely  to  survive  as  it  ages,  a multivariate  generalization  of  this  is  that 


IUjLX l 

Fix) 


decreases  in  x > 0 for  all  t > 0. 


(3.1) 


This  condition,  however,  is  ambiguous  since  x can  increase  in  various  ways  and 
^ can  he  restricted  in  certain  ways.  Several  versions  of  (3.1)  will  now  he 
given.  It  is  also  assumed  that  similar  conditions  hold  for  all  marginal 
distributions.  Here  1 **  (1 1). 


(a)  lAl  decreases  in  x ' 0 for  all  t > 0. 

V(Vl) 


(b) 


F& L+ JJL 

F(xl) 


decreases  in  x > 0 for  all  t s 0. 


(c) 


F(a+,  tj) 

F(x) 


decreases  in  x > 0 for  all  t N 0. 


(d)  decreases  in  x > 0 for  all  t > 0. 

V(x) 
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Another  version  in  a similar  spirit,  although  not  exactly  in  the  same  form 
as  (3.1)  is 

F((xt  + t)t,  x) 

(e)  for  i - 1,...,  n,  — decreases  in  x.  > 0 for  all  t > 0 

F(x) 

and  all  >_  0,  j i . 

Here  the  notation  ((x  + t)  , x)  means  (x  , , . . , x , , x + t,  x x ). 

i i ~ i i-l  i i+1  n 

Version  (c)  is  due  to  Harris  [17]  and  to  Brindley  and  Thompson  [11];  (a),  (b) 
and  (d)  are  due  to  Marshall  [21];  and  (e)  (in  another  form)  is  due  to  Johnson 
and  Kotz  [20].  Many  other  versions  are  possible,  but  it  appears  that  of  all  of 
these  (c)  is  perhaps  the  most  important.  One  reason  for  this  is  that  (c) 
captures  the  idea  of  the  original  model,  in  that  time  for  the  components  runs 
at  the  same  rate  (i.e.  t = t^  = •••  = tR) . Furthermore  the  different  x^,...,  x^ 
allows  for  the  possibility  that  the  components  are  of  different  ages,  as  is 
often  the  case  in  practice.  Concept  (c)  is  designated  MIFR  and  is  entensively 
discussed  in  Barlow  and  Proschan  [3].  A dual  concept  with  increasing  replacing 
decreasing  is  called  MDFR.  Other  reasons  for  the  importance  of  MIFR  are  the 
following  properties. 

(i)  A univariate  MIFR  lifetime  is  IFR. 

(ii)  The  union  of  independent  sets  of  MIFR  lifetimes  is  MIFR. 

(iii)  The  marginal  lifetimes  of  MIFR  lifetimes  are  MIFR. 

(iv)  Series  systems  formed  with  MIFR  lifetimes  are  MIFR. 

(v)  A lifetime  is  MIFR  and  MDFR  if  and  only  if  it  has  the  multivariate 
exponential  distribution  known  as  the  MVE  (see  Chapter  5 of  Barlow  and 
Proschan  [3]). 

The  second  notion  upon  which  other  versions  of  multivariate  IFR  has  been 
based  is  a generalized  form  of  (2.4).  These  versions  were  given  by  Marshall  [21] 
and  involve  the  idea  that  log  F (x  , . . . , xr)  is  concave  in  some  sense.  This 


i 


concavitv  can  be  taken  in  the  usual  sense  or  along  various  curves  ami  linos  in 


Rn.  The  properties  of  the  various  classes  which  emerge  have  been  briefly 
discussed  by  Marshall  £j2 lj  . 


4.  Multivariate  I FRA . The  central  theoretical  role  played  bv  the?  univariate 
TFRA  class  derives  from  the  properties  satisfied  hv  thin  class  (see  i)  - iv) 
following  (2.7))  rather  than  from  the  fact  that  it  has  an  increasing  average 
failure  rate.  It  is  not  surprising  then,  that  various  multivariate  versions  of 
IFRA  have  been  defined  using  a generalization  of  one  of  these  basic  properties. 
The  classes  of  lifetimes  which  result  from  these  generalizations,  however,  all 
fail  to  satisfy  some  basic  property.  Block  and  Savit.s  [8]  have  proposed  a class 
of  multivariate  TFRA  lifetimes  which  is  based  on  a mathematical  property  of 
univariate  IFRA  lifetimes.  This  class  contains  a rich  variety  of  multivariate 
lifetimes  and  also  satisfies  all  of  the  fundamental  properties  which  one  would 
expect  for  a multivariate  IFRA  class. 

The  class  of  multivariate  lifetimes  proposed  by  Block  and  Savits  is  giver 
in  the  following  definition  which  is  a generalization  of  the  univariate  property 
(2.7). 


Definition  4.1. 


random  vector  T 


T.et  T =•  (T.j , . . . , T ) be  a nonnegative 
is  said  to  be  MI FRA  if 

Ea  [h(T)l  £ F.[ha(T/«)1 


random  lifetime. 


The 


(4.1) 


for  all  continuous  nonnegative  nondecreasing  functionsh  and  all  0 < a < 1. 

Other  conditions  which  have  been  proposed  for  multivariate  TFRA  are  given 
in  the  next  definition.  As  in  the  previous  material,  we  rely  heavily  on  the 
notation  and  terminology  of  Barlow  and  Proschan  [7)  with  the  one  oxccpHon 
mentioned  in  Section  1.  The  life  function  t corresponding  to  >i  ;-tr-irtu»e  function 
$ is  called  monotone  (coherent)  if  if  is  monotone  (coherenl ) . See  Esary  and 
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Marshall  [15]  for  a discussion  of  life  functions. 

Definition  4.2.  Let  T ■ (T T ) be  a nonnegative  random  lifetime  with 

—in 

survival  function  F(t)  - P(T  >>  t)  where  ">>"  means  strict  inequality  holds 

for  each  component.  The  vector  T is  said  to  satisfy  condition  if  the 

condition  following  letter  is  satisfied: 

A:  F01  ( t ) <_  F(at)  for  all  0 < a £ 1 and  all  0 <_  _t. 

B:  T is  such  that  each  monotone  system  formed  from  T is  univariate  IFRA. 

C:  T is  such  that  there  exist  independent  IFRA  random  variables  X^,...,  X^ 

and  monotone  life  functions  t^,  i ° 1,...,  m such  that  (X^,...,  X^) 

for  i *■  1, . . . , m. 


T is  such  that  there  exist  independent  IFRA  random  variables  X^,...,  X 

and  nonempty  sets  S.  of  {1,...,  k)  such  that  T.  £ X.  for  i * 1,...,  m. 

jeSi  3 

D:  T is  such  that  there  exist  independent  IFRA  random  variables  Xj^...,  X^ 

and  nonempty  subsets  S.  of  {1,...,  k)  such  that  T.  ■*  min  X for  i * 1,...,  m. 

jeS.  3 

E:  T is  such  that  the  minimum  of  any  subfamily  of  T T is  IFRA. 

— — 1 m 

F:  T is  such  that  min  a T is  IFRA  for  all  a > 0,  i ■ 1,...,  m. 

i 1 1 1 ” 

Conditions  A,B,C,D,E,F  have  been  given  by  Esary  and  Marshall  [16]  and  condition  £ 
was  given  by  Block  and  Savits  [9]. 

The  following  relationships  hold  between  MIFRA  and  the  seven  conditions 
(see  Block  and  Savits  [10]  for  proofs  and  examples). 


Q 


» l 

U // 

\ / 

Ml  FRA 

A(-)F  R 


Figure  1 

With  the  exception  of  the  implication  ) C (which  is  conjectured  to  he  false). 
Block  and  Savlts  (10)  have  shown  that  no  other  implications  are  possible. 

The  following  properties  are  basic  for  a class  of  multivariate  1FRA 
distributions . 

(PI):  Closure  under  the  formation  of  monotone  systems. 

(P2) : Closure  under  limits  in  distribution. 

(P3) : Marginals  are  in  the  same  class. 

(P4) : Closure  under con junction  of  independent  lifetimes. 

(P5) : Closure  under  scaling. 

(P6) : Closure  under  convolution. 

It  has  been  shown  by  Block  and  Savits  [81  that  the  MIFRA  class  satisfies  all 
six  of  these.  Furthermore  in  a subsequent  paper,  Block  and  Savits  [101  have 
demonstrated  that  each  of  the  seven  conditions  in  Definition  4.2  fail  to  satisfy 
at  least  one  of  these  properties.  Specifically  A,£,  E and  F fail  to  satisfy 
(PI);  B and  C fail  to  satisfy  (T5);  and  D fails  to  satisfy  (Pb)  . 

5.  Multivariate  NBU  and  NBUE 

Unlike  the  previous  two  multivariate  classes  discussed,  only  preliminary 
work  has  been  done  in  the  situation  of  multivariate  NBU  and  NBUE  classes.  We 
shall  present  12  possible  definitions  of  NBU  and  a similar  group  for  NRUE. 
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The  only  versions  of  multivariate  NBU  of  which  we  are  aware  in  the 
literature  have  been  discussed  by  Buchanan  and  Singpurwalla  [121  and  Marshall 
and  Shaked  [221.  These  generalize  (2.7),  have  a similar  interpretation,  and 
are  related  to  condition  (3.1).  They  specify  that  for  the  random  lifetime 

i - <T1 v- 

p{T  >>  x + Jt  | T >>  xl  “ ^ 1 F(t)  - P(T  » _t}  (5.1) 

for  various  choices  of  x and  Conditions  1)  - 3)  in  the  following  are  the 
ones  given  by  Buchanan  and  Singpurwalla.  Condition  4)  is  of  the  same  type  and 
was  given  by  Marshall  and  Shaked  [22]  and  Condition  5)  is  similar  in  spirit  to 
e)  of  Section  3. 


1)  F((x  + t)l)  < F(xl)  F(tl)  for  all  x > 0,  t > 0. 

2)  F(xl.  + t)  < F(xl)  F(t)  for  all  x > 0,  t > 0. 

3)  F(x  + t)  < F(x)  F(t)  for  all  x>  0,  1 1 0. 

4)  F(x  + t)  <_  F(x)  F(t)  for  all  x > 0,  t > 0 which  satisfy  (x^  - x^)(t^  - t^)  0 

for  all  i,  j. 

5)  F(x  + tEj)  < F(x)  F(e^)  for  all  x > 0 and  i = 1,...,  n where 

■ (0,...,  0,  1,  0,...,  0)  and  1 appears  in  the  ith  position. 


The  next  condition  is  a variant  of  these  and  actually  specifies  several  conditions 
depending  on  the  class  of  life  functions  chosen. 

6)  P{t(T  - x)  > t |T  > x}  <_  P{r(T)  > t}  for  all  t 0,  x > 0 and  all  t in 
a class  of  life  functions  (e.g.  for  all  minimums). 


If  the  class  of  life  functions  specified  is  the  class  of  minimums  (i.e.  series 
systems)  then  6)  can  be  interpreted  as  the  probability  of  survival  for  more  than 
t units  of  any  series  system  formed  from  a set  of  components  of  age  at  least  x 
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being  less  than  the  probability  that  the  corresponding  new  series  system 
survives  more  than  t units. 

Conditions  7)  - 9),  to  be  given,  are  similar  in  spirit  to  the  conditions 
B)  - F)  in  Section  4. 

7)  t(T  , T^)  NBU  for  all  t in  a certain  class  of  life  functions  (e.g.  for 
all  monotone  functions,  all  minlmums,  all  sums). 

8)  There  exist  independent  NBU  Xj.,...,  Xfc  and  monotone  life  functions 

i = 1 n such  that  T = x1(X)  for  i = 1,...,  n. 

9)  x(a^T^)  NBU  for  all  a^  0 and  all  t in  a certain  class  of  life  functions. 

Again  each  of  7)  - 9)  represent  several  conditions  depending  on  the  specification 
of  the  or  t. 

The  following  condition  is  the  generalization  of  a characterization  of  NBU 
given  in  Block  and  Savits  [7]. 

10)  E(h(T  - x)  | T > x)  £ E(h(T))  for  all  nonnegative,  nondecreasing  and 
continuous  h. 

Another  possible  generalization  of  NBU  is  through  the  concept  of  multivariate 
shock  models  of  a type  discussed  in  the  univariate  case  by  Barlow  and  Proschan 
13].  Several  multivariate  shock  models  have  been  discussed  by  Marshall  and 
Shaked  [22]  and  by  Block  [4].  Using  the  notation  of  Section  3.0  of  the  latter 
paper  we  have 

N N 

11)  T *=  ( 1 Xjj...,  J Xin)  where  (X^,...,  X^)  has  independent  exponential 

marginals  and  these  vectors  are  independent  for  i ■=  1,2,...  and  N is 
univariate  NBU  and  independent  of  the  X^. 

There  are  many  variants  of  this. 


A final  version  of  multivariate  NBU  can  be  given  by  generalising  a 
characterisation  of  NBU  given  by  Block  and  Savits  (6]  and  using  the  Laplace 
transform.  This  involves  giving  a discrete  condition  analogous  to  1)  - 5) 
on  coefficients  related  to  the  Laplace  transform  of  a random  lifetime  T.  Since 
this  is  notationally  involved  we  will  omit  details. 

Comments  on  various  possible  NBUE  definitions  can  similarly  be  made.  Various 
Integrated  versions  of  1)  - 6)  can  be  given.  Four  versions  of  1)  - 3)  are 
given  by  Buchanan  and  Singpurwalla  [12].  Versions  of  7)  - 9)  can  be  given  where 
the  t and  are  specifically  taken  to  be  sums.  Since  a characterisation  of 
NBUE  is  given  in  Block  and  Savits  [7]  a multivariate  version  analogous  to  10)  con 
be  given.  If  in  11)  N is  assumed  to  be  NBUE  instead  of  NBU  a multivariate  NBUE 
definition  is  obtained.  Furthermore  the  comments  in  the  paragraph  following  11) 
apply  as  well  to  the  NBUE  case. 

As  mentioned  previously,  little  work  has  been  done  on  multivariate  NBU  and 
NBUE  classes.  Determining  which  definitions  are  most  fundamental  and  how  these 
concepts  are  related  to  one  another  still  remains  to  be  done.  This  problem  is 
currently  being  studied  by  the  authors. 
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